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Abstract 

In this paper we consider a bistable reaction-diffusion equation in 
unbounded domains and we investigate the existence of propagation 
phenomena, possibly partial, in some direction or, on the contrary, 
of blocking phenomena. We start by proving the well-posedness of 
the problem. Then we prove that when the domain has a decreasing 
cross section with respect to the direction of propagation there is com¬ 
plete propagation. Further, we prove that the wave can be blocked 
as it comes to an abrupt geometry change. Finally we discuss various 
general geometrical properties that ensure either partial or complete 
invasion by 1. In particular, we show that in a domain that is “star- 
shaped” with respect to an axis, there is complete invasion by 1. 
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Introduction 


1.1 The problem 

In this paper we consider the following parabolic problem 

jd t u(t, x) — Au(t, x) = f(u(t, x)), for t E M, x G 12, 

\^d v u(t,x) = 0, for t e M, x G <912, 

where 12 is an unbounded domain such that 

12 = {(xi, x'), X\ G M, x' £ ca(xi) C . (1.2) 

Throughout the paper, we assume to not to depend on x\ for x\ < 0 and 
/ is a bistable nonlinearity. We call such domains “cylinder-like”. We give 
more precisions on the general assumptions in the next section but one can 
see in Figure [l] various examples of the type of domains 12 we have in mind 
in dimension 2. 



Figure 1: Examples of “cylinder-like” domains 12. Propagation is considered 
in the direction of increasing X\. 


The aim of this paper is to describe invasions coming from — oo. That is, 
waves propagating from x\ = — oo in the direction of increasing X\. These 
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waves have implications for the initial value problem for initial data with, 
say, compact support. We can also envision domains with two different axes, 
with the direction of propagation changing as t —> —oo and as t —> + 00 . 
Think for instance of “bent cylinders”. But for simplicity, in this paper we 
only present the results in the framework of a single axis. The main results 
in this paper are the following. 


We first study the existence of transition fronts in this setting, i.e the 
existence and uniqueness of an entire solution u of the parabolic prob¬ 
lem (1.1), such that 

u(t, x ) — p(x 1 


cl) —y 0 as t —y —00 uniformly in fl, 


(1.3) 


where ( ip , c) is the unique bistable travelling front solution of: 

f p” + op' + f(<p) = 0 in M, 

1y9(-oo) = 1, <^(+ 00 ) = 0, </?(0) = 9. 


• We establish the blocking of the solution if there exists a G R such that 
the measure of 11 fl {a < X\ < a + 1} is small, 

• We prove that there is axial partial propagation when 11 contains a 
straight cylinder in the aq-direction, whose cross section is a ball of 
sufficiently large radius. 

• We give conditions on the geometry of hi under which there is complete 
invasion that is, u(t, x) —* 1 as t —* 00 at every point x: 

■ In the case of general domains with decreasing cross section, 

• Under the condition of partial convexity and domains with in¬ 
creasing cross section, 

■ For domains “star-shaped” with respect to an axis. 


Let us specify some definitions on propagation and blocking properties. Let 
u 00 (x) := lim u(t,x ) for all (see Theorem 1.4 in the next section). 

t —>-+oc 


Definition 1.1 (Blocking and partial propagation). - Let u be the solution 
of problem (1.1) satisfying (1.3). We say that u is blocked when 


Uoo(x i,x') —>■ 0 as X\ —> + 00 , (xi,x') G U. 


Otherwise we say that there is partial propagation in fl. 
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Definition 1.2 (Axial partial propagation). - Let u be the solution of prob¬ 
lem (1.1) satisfying (1.3), we say that there is axial partial propagation if 


inf Uoo > 0, 

RxBjj 


for some R > 0 and where B' R is the ball of radius R centred at 0 in 1 . 


Definition 1.3 (Complete invasion). - Let u be the solution of problem (1.1) 
satisfying (1.3), we say that there is complete invasion if = 1 in f l. 


This problem is of interest in different domains of application. In pop¬ 
ulation dynamics, to study the evolution of a population going through an 
isthmus or for fish populations going through a straight for example. In med¬ 
ical sciences to model the evolution of a depolarisation wave in the brain or to 
study ventricular fibrillations. These motivations from biology and medicine 
are detailed in section [2j where we also discuss the interpretations of our 
results in these contexts. 


1.2 Previous mathematical results 

In this paper, we are interested in the existence and properties of invasion 
fronts in the direction X\ (see [5] for definition) for a bistable reaction- 
diffusion equation in unbounded domains. This problem was partly studied 
in m in special cases. There the authors proved that in two or three di¬ 
mensions, if Vi is the succession of two rectangles (or parallelepipeds) with 
different widths r and R , then one can find conditions on r and R such that 
the solution is blocked when it goes from the small rectangle to the large 
rectangle. The authors used a specific symmetrization method that does not 
work for other domains. In [7], the authors proved the existence of gener¬ 
alised transition fronts in exterior domains, which pass an obstacle. They 
analysed the interactions between travelling wave solutions and different kind 
of obstacles. In particular, they gave geometrical conditions on the shapes 
of the obstacles under which there is complete invasion by the most stable 
state. They also constructed an obstacle that blocks the bistable travelling 
wave in some area of the domain. In this case, there is only partial invasion. 

In this paper, we investigate the question of propagation/blocking prop¬ 
erties for rather general “cylinder-like” geometries. Roughly speaking, by 
this we mean to say that there is an axis, say the aq-axis, in the direction 
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of which the domain is unbounded (both as ii —* — oo and x\ —> +oo). 
Furthermore, for simplicity, we assume that the domain is a genuine cylinder 
in the part {aq < 0}. We first construct a generalised transition front using 
the techniques of [7J. Then the main purpose of this paper is to analyse the 
geometrical conditions, under which, the wave is blocked or propagates and, 
in the latter case, whether the invasion is partial or complete. 

Related questions have been studied previously in the literature. In EU 
Grindrod and Lewis approximate the normal speed of the interface of a front 
using the Eikonal approximation and numerical simulations to point out that 
the presence of a symmetric narrowing and widening will decrease the av¬ 
erage speed and that an abrupt widening could lead this speed to 0. More 
recently in [30!, 28j, Lou, Matano and Nakamura study the effect of undu¬ 
lating boundary on the normal speed of the interface of a front analysing a 
curvature-driven motion of plane curves in dimension 2. They give necessary 
and sufficient conditions on the undulations of the domain for the existence 
of generalised travelling waves in this context. 

1.3 Notations and assumptions 

Here are some notations that we will use throughout this paper. First M - = 
(—oo,0) and M + = (0,+oo). The space dimension is N > 2 and we write 
the space variable as x = (aq, x') G R N with aq G M and x' G M^ -1 . 

We will use balls both in as in Thus to distinguish the two 

cases, we use the notation 

Br(x o) := {a: G R N , \x — Xo| < -R} 

for the ball of radius R and centred at xq in whereas 

B' r {x o) := {x' G R n ~\ \x' - z 0 | < R } 

denotes the ball of radius R and centred at xq in M^ -1 . If xq = 0, we will 
simply write Br and B' R . The outward unit normal at x G dfl will be denoted 

u(x) = (vi(x), v\x)) Glx M 7 ^ 1 . 

Now let us state precisely the conditions on the bistable nonlinearity /. First 
we assume that / G C 1,:L ([0,1]) and there exists 6 G (0,1) such that 

/(0) = m = /(1) = 0, /'(0) < 0, /'(1) < 0, (1.5a) 
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f(s) < 0 for all s G (0,0), f(s ) > 0 for all s G (9, 1). (1.5b) 

Moreover we assume that 1 is “more stable” than 0 in the sense that: 


f(s)ds > 0. 


(1.6) 


We let (3 G (9 ,1) denotes the real number: 

f3 = inf G [0,1], J f(s)ds > 0 j- . (1.7) 

As mentioned before, (ip, c ) will always denote the one-dimension travelling 


front with <£>(0) = 9, i.e. (<p,c) is solution of (1.4). 

Concerning the domain 12 C , we have in mind to consider cylinders 
around the ap-axis with varying cross section but the various results are valid 
for much more general domains and our only general assumptions are that 


12 is a uniformly C 2,a domain of M JV and 


N 


( 1 . 8 ) 


Q n { 


x G 


D N 


, xi < 0} = M Xu where u> is an open subset of 1 N 1 . (1.9) 


Condition (1.9) means that the portion of lying in a half space (here 
{xi < 0}) is a straight cylinder. Thus, the change in the geometry of the 
domain takes place in the half space M + x M 7 ^ 1 . 

1.4 Main results 

We first prove the wcll-posedness of the problem with the following theorem. 


Theorem 1.4. Let 12 C 


t,N 


and f satisfy conditions (1.8)-(1.9) and (1.5) 


(1.6), let (<p,c) be the unique solution of (1.4), with c > 0. Then there exists 


a unique entire solution of 


dtu(t,x) — Au(t,x) = f(u(t,x)), for t G M, x G 12, 
d v u(t, x) = 0, for t G M, x G <912, 


such that 


u(t, x) — <p(x i — ct) —> 0 as t —>■ —oo uniformly in 12. 
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Moreover, u satisfies u t (t,x ) > 0, 0 < u(t,x) < 1 for all (i,i)eKxU and 
lim u(t ,.) = Moo m Cf oc (Ul) where Uoo is stationary solution, that is, verifies 

t—t+oa 


AMoo + /(Moo) =0 X G ft, 
dv^oo = 0 x G <912. 


( 1 . 10 ) 


This theorem is essentially due to {7j, who considered an exterior domain. 
But their proof readily extends to our case. We sketch it in section 3. 
Remark 1.5. This result also applies to bent and right-ended cylinders as 
those illustrated in section 13.41 

Most of this paper is devoted to analyse the large time behaviour of u(t, x), 
that is understanding when it is blocked, when there is partial propagation or 
axial partial propagation or when there is complete invasion. We start with 
the case of domains with decreasing cross section for which we can prove that 
there is complete invasion. 


Theorem 1.6 (Complete propagation into decreasing domains). Assume 
that for all x G <912, isi(x) > 0, where vfix) is the first component of the 
outward unit normal at x. Under the same assumptions on 12, / and (yyc) 
as in Theorem l.f, the solution u of (1.1) satisfying (1.3) propagates to 1 in 
12 in the sense of Definition \1.3\ That is 


u(t, ■) —> 1 as t —> +oo in 12. 


Moreover, if we assume that 

12 fl {x G M. N , X\ > 2} = (l, Too) x ay, 

for some l > 0 , ay C M^ -1 , then c is the asymptotic speed of propagation, i.e 


V c > c, lim sup u(t,x) = 0, 

t->+oc Xl> gt 

V c < c, lim inf u(t,x) = 1. 

t —>- + OC xi <ct 

This theorem means that there exists a generalised transition front rep¬ 
resenting an invasion of 0 by 1 for problem (1.1). 

Remark 1.7. If we assume that 12 is a cylinder-like domain with axis on the x\ 
axis (i.e w(-) defined in (1.2) is open, connected, smooth and bounded), the 
condition on vfix) is fulfilled as soon as the cross section is non-increasing 
with respect to xi, but it applies to more general domains (see the upper 
middle figure of Figure [l] for an example in dimension 2). 











Next we study the existence of blocking phenomena and prove the follow¬ 


ing 

Theorem 1.8 (Blocking by a narrow passage). Under the same assumptions 
on 12, / and (<p,c) as in Theorem l.f and let a,b be fixed with —00 < a < 
b < Too. There exists e > 0 small enough depending on 

12 fl {x G R n , b < x\ < b T 1} (1.11) 

such that if 

112 fl (x G M. N , Xi G (a, 6)} | < e, (1.12) 

then the unique solution u of fll.l ) satisfying (1.3) is blocked in the right part 
of the domain in the sense of Definition 1.1, That is 


u(t, •) -r Moo in 12 as t -T Too with u^x) —* 0 as x\ —* Too. 


A particular case of this theorem is the example of abruptly opening 
cylinders in the direction of propagation x\. More precisely, for a fixed right 
hand half of the domain 12 + := 12fl{xi > 0}, there exists e > 0 small enough 
such that if 12fl{xi < 0} C M x B e , the advancing wave is blocked. From the 
proof, essentially, we will see that it is blocked at the abrupt passage from 


establishes that for a large left cylinder, a generalised travelling front can 
never be blocked. Thus rescaling the space variable to fix a given measure 
to the left cylinder and enlarge the right side of the domain does not neces¬ 
sarily work. Figure [2] displays examples of 2D domains where this blocking 
phenomenon arises. 

Remark 1.9. Notice that rescaling / into A/ for some A positive, then e 
becomes —^=, which means that if we increase the amplitude of / then, 

the threshold £ decreases and blocking will take place in thinner channels. 
Furthermore, if we add a diffusion coefficient D > 0 in front of the second 
order term (scaled to 1 in this paper) then e becomes \ [De, which means 
that if D increases then the solution is blocked for larger diameters. 

The blocking phenomenon as well as these scalings can be interpreted in 
the following manner. As the species moves through a narrow passage, it 
all of a sudden emerges into a widen open space. As it emerges, diffusion 
disperses the species, possibly to lower values of density where the reaction 
becomes adverse. Therefore, too narrow a passage does not allow for the 


a narrow cylinder into a wide domain. On the other hand, Theorem 1.13 
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Figure 2: Examples of domains in which the wave is blocked at the passage 
into the wide region, for small enough e > 0. 


species to reconstitute a strong enough basis to invade the right hand half of 
the domain and the invasion wave is stuck at this passage. Clearly, a large 
diffusion accentuates this problem, while a stronger /, by allowing large 
densities to transit through the channel has the opposite effect. 

Now considering domains that are not narrow, we show that this does not 
happen if the cross section is large enough. We prove that when 12 contains 
a straight cylinder in the ap-direction of radius large enough, there always is 
propagation. 


Theorem 1.10 (Axial partial propagation). Under the same assumptions as 
in Theorem 1 . 4 , there exists Rq > 0 such that i/lx B' R C 12 for an R > Rq, 
the unique solution u of (1.1) satisfying (1.3), propagates along the x\-axis 
in the sense of Definition \1.2\ That is 


u(t , •) —x Woo in 12 as t —x +00 with inf > 0. 

RxB' r 

We give more details about how this constant is determined in section |6j 

Remark 1.11. Notice that in this theorem the invasion may be partial. Indeed 
depending on the shape of 12, there could be subdomains where the solution 
Moo is close to 0 (see [TUTO] and Figure[9]in section[6]of this paper for examples 
of domains where it happens). 

I 11 the following theorems, we focus on the geometry 12 such that the 
invasion is complete, i.e. Uoa = 1. The next theorems are proved using a 
sliding method (either sliding a ball or a cylinder), which is inspired from |[9j. 
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Theorem 1.12 (Complete propagation in “star-shaped” domains). We as¬ 
sume that the conditions of Theorem \1.10\ are fulfilled. Moreover we assume 
that at each point on the boundary x = {x\,x') G <912, the outward unit normal 
v makes a non-negative angle with the direction x'. More precisely, writing 
v = (z/ 1; u') (with v' G we assume: 

u' ■ x' > 0 for all points x = (x\,x') G <912. (1-13) 


Then, the invasion is complete, namely, u 00 = 1 in 12. 


We call a domain that satisfies the property (1.13) a star-shaped domain 
with respect to the x\-axis. One can look at Figure [l| (upper left, lower left 
and middle) for examples of such domains in dimension 2. More examples 
and counterexamples of such domains are given at the end of section [7j 
The previous Theorem covers a large range of geometries for 0, however 
the following Theorem completes this result for another type of geometry. It 
relies on a different proof. 


Theorem 1.13 (Complete propagation into increasing domains). Under the 
same assumptions as in Theorem\TA there exists R\ > 0 such that if 


M x B' r C 12, for R > Ri, 

Q r L := {(x\,x') G 12, X\ > L} is convex, 

{x G 12, Xi < L + R} C M. x Be, 

Vx = (xi, x') G <912 x\ < L + R => vfix) < 0 


for some L,C > 0, then the invasion is complete in the sense of Definition 


unit normal at x. 


1.3. namely Uoo = 1 in 12. Here ufix) is the first component of the outward 


In the proof of this theorem we use a solution of the stationary elliptic 
equation in a ball with Dirichlet conditions. We then slide this ball in (M x 
B' r ) U 12^ to get a lower estimate on the solution u ^ in this domain and 
then conclude with the maximum principle and Hopf Lemma. This proof 
could be adapted to somewhat more general domains as soon as it satisfies a 
sliding ball assumption, i.e every point in the domain could be touch, from 
within the domain, by a ball of radius R — 6 for some 6 > 0 such that for 
x G <912 fl B r , u(x) ■ x > 0. 

We thus point out two different types of behaviour of the solution in the 
case of widening cylinders. If the diameter of the cylinder before the change 
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in geometry is small and the change is abrupt, we prove that there is blocking, 
whereas if the cylinder is wide enough before the change of geometry there 
is propagation even if the domain abruptly increases. One question that 
remains open is the existence of a sharp threshold, when considering, say 
a monotone 1-parameter family of domains. The question is thus to know 
whether the solution is blocked below a critical threshold and propagates 
above it. 

In this paper we consider mainly generalised transition waves. But of 
course our results and methods bear important consequences for associated 
Cauchy problems. 

Corollary 1.14 (The associated Cauchy problem). Let u be the solution of 
the following Cauchy problem 

{ d t u(t, x) — Au(t, x) = f(u(t, x)), for t G (0, +oo), x G 12, 

d v u(t , x) = 0, for t G (0, Too, ) x G <912, 

u(0, x) — u 0 (x), for x E ft. 


Then 


Under the same assumptions as in Theorem 1.8 on ft, f and (</?, c), the 
propagation of the solution is blocked for any u 0 , of which the support is 
contained in the left part of the domain , i.e supp{u 0 } C (x G 12, X\ < 
a}, 


■ Under the same assumptions as in either Theorem m Theorem \1.10\ 
Theorem 1.12 or Theorem 1.13 on 12, f and (<p,c), the propagation 
results would still hold if the domain {x G 12, u 0 (x ) >6 + 5} is large 
enough, for some 5 > 0. 

This article is organised as follows. In section [2] we detail the motivations 
from biology and medicine behind this problem. We indicate the interpreta¬ 
tions of our hirelings in these contexts. We prove the existence of transition 
fronts (Theorem 1.4) in section [3j We devote section [4] to the proof of The¬ 
orem |1.6 on the propagation of fronts in decreasing domains. Section [5] is 


devoted to the proof of Theorem |1.8[ that is the blocking of the front if 12 
contains a narrow passage. Then we derive a sufficient condition under which 
there is (possibly partial) propagation in section [bj There, we prove Tlieo- 
Next we turn to further complete propagation results. We prove 


rem 


1.10 


Theorem 11.121 in section [7] and Theorem 11.131 in section |H1 


12 


















2 Motivations from medicine and biology and 
interpretations of the results 


Equation (1.1) is classical in physics, for the representation of phase transi¬ 
tions, and in biology. A particularly important motivation for studying the 
effect of geometry on propagation of waves in equation (1.1) comes from the 
modelling of Cortical Spreading Depressions (CSDs) in the brain. These are 
transient and large depolarisations of the membrane of neurons that propa¬ 
gate slowly (with a speed of the order of 3mm/min) in the brain. CSDs are 
due to abnormal ionic exchanges between the intra- and extra-cellular space 
of the neuronal body that slowly diffuse in the brain. There are two stable 
states, the normal polarised rest state and the totally depolarised state, with 
a threshold on the ionic disturbances for passage from one stable state to the 
other. After a depolarisation lasting 3 or 4 minutes, several mechanisms take 
place that repolarise the neurons. 

Mathematical studies of CSDs focus generally on the depolarisation phase. 
Having in view this phase, CSDs are modelled by a bistable reaction-diffusion 
equation like (1.1), where 0 < u < 1 represents the degree of depolarisation. 
In this context, u = 0 stands for the normal polarised state and u = 1 repre¬ 
sents the completely depolarised state. Thus, a depolarisation is represented 
by an invasion of the state u = 1. The complete depolarisation mechanism 
only takes place in the grey matter - the part of the brain that contains the 
neuronal bodies - while simple ionic disturbances diffuse everywhere and are 
absorbed in the white matter - the part of the brain where no neuron body 
and only axons are to be found. In this context, the domain fl in equation 
(1.1) can be thought of as representing a portion of grey matter of the brain. 
The Neumann boundary conditions here is a classical simplification resulting 
from considering that the grey matter is isolated. 

The question of existence of and mechanisms to break CSDs in human 
brain is of great importance for the understanding of strokes or migraines 
with aura among other pathologies. CSDs were first observed by Leao in 1944 
m When neurons are depolarised, they cannot propagate nerve impulses 
and this generates several symptoms. For example, CSDs are suspected of 
being responsible for the aura during migraines with aura [2411251 [22] . The 
aura is a set of hallucinations, mainly visual and sensory, that appears during 
some migraine attacks. In ischemic strokes in rodents, it has been proved 
that each CSD increases the neurological damage by approximately 30% [33] 
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and therapies aiming at blocking the appearance of CSDs have shown very 
promising results in rodents [161 [35]. Unfortunately they turned out to be 
inefficient in humans. 

Moreover there are obvious difficulties in observing CSDs in the human 
brain. These difficulties are related to the complex geometry of the human 
cortex. With EEG, each electrode records brain electrical activity from a cor¬ 
tical area of several square centimeters whereas the depolarized area during 
CSDs is much smaller (around few square milimeters). Alternative medical 
imaging methods record accompagnying symptoms of CSDs (cerebral blood 
flow increase or decrease) but not the depolarization itself. CSDs can be 
recorded only in acutely injured human brain undergoing neurosurgery, using 
cortical electrode strips similar to the technology applied in epilepsy surgery. 
But in such a pathological situation, the brain reacts quite differently from 
a healthy brain. 

For many years, it was believed that CSDs were an artifact produced in 
animal experiments and without significance for human neurological con¬ 
ditions. Their occurence in human cortex was thus a matter of debate 
[H H91 EH [321 Sol [39113]. However in the recent years, a medical consensus 
considering that CSDs occur in human brain has gained ground (see [26] and 
references therein). Consequently the question of how to block these CSDs 
became important. For migraines with aura, CSDs have been observed to 
become extinct when propagating into a sulcus HU, that is in the bottom 
of a convolution. Generally the grey matter is thicker in a sulcus than in 
the rest of the cortex. For a CSD during a stroke, understanding how far 
the CSD will propagate in a specific patient could help in anticipating the 
dammage caused by the stroke, ft could also help in the selection of patients 
for future clinical trial. 

Mathematical modelling may help in understanding these phenomena 
since experiments on rodents cannot be conclusive for humans due to the 
differences of the brain morphologies. The grey matter of the human brain 
composes a thin layer at the periphery of the brain with its thickness subject 
to abrupt and large variations. The rest of the human brain is composed 
of white matter. On the opposite, the rodent brain is almost round and 
entirely composed of grey matter (see Figure [3j) . The variations of the grey 
matter thickness in the human brain could explain the extinction of CSDs. 
Indeed, as we establish it here, very abrupt variations of the domain may 
block propagation 

Chapuisat and Grenier D33 first stated the hypothesis that the morphol- 
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Figure 3: Morphology of the brain, cross sections. Left, human brain: grey 
matter is only a thin layer at the periphery. Right, rodent brain: it is almost 
entirely composed of grey matter. 

These images are from the University of Wisconsin and Michigan State 
Comparative Mammalian Brain Collections, and from the National 

at the following website 
Preparation of all these 
images and specimens have been funded by the National Science Foundation, as 
well as by the National Institutes of Health. 


Museum of Health and Medicine, available 
http: / / www. br ainmuseum. or g / sections / index, ht ml, 


ogy of the human brain may prevent the propagation of CSDs over large 
distances, unifying several observations by biologists. Several effects may be 
combined. First, variations of the grey matter thickness may prevent prop¬ 
agation of CSDs. This point is studied numerically in m and proved in a 
special case in H3j. Furthermore, the absorbing effect of the white matter 
may stop the propagation of CSDs in areas of the brain where the grey mat¬ 
ter is too thin. This phenomenon is studied in [13] and more recently in [[TJ. 
Finally the convolution of the human brain may also enhance the absorption 
effect of the white matter and thus the blocking of the CSDs. This point is 
studied numerically in [S6j. Here we prove in a general setting that abrupt 
transitions from thin to large regions will block waves. 

These results can be interpreted in a biological framework. If we assume 
that the mechanisms triggering CSDs are the same in humans and in rodents, 
CSDs could be initiated in the human brain but only in areas where the grey 
matter is large enough and these CSDs will not propagate over large distances 
since they will be blocked by a sudden enlargement of the grey matter, or 
if the grey matter becomes too thin as compared to the white matter, or if 
they go through a sharp convolution of the brain. On the opposite, in rodent 
brains that are like balls of grey matter, CSDs could be triggered anywhere. 
This may explain why it has been difficult to observe CSDs in human brains 
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during a stroke. 

In the case of migraines with aura, the symptoms of aura correspond to a 
dysfunction of certain neurons. When translating the symptoms described by 
patients on a map of the cerebral areas, it is observed that the dysfunction 
of neurons propagate at a speed of 3mm/min which makes biologists and 
physicians think that the aura is due to CSDs. It has been shown that the 
aura always stops at the bottom of a very deep sulkus. This phenomenon 
has been studied numerically in [20] and it strengthens the hypothesis of 
the role of the morphology of the human brain in the propagation of CSDs. 
However a complete theoretical study of the influence of the geometry of the 
grey matter on the propagation of CSDs has never been carried out up to 
now. 

Another motivation for the study of the effect of geomery on propaga¬ 
tion in equation (1.1) comes from the study of ventricular fibrillations where 
blocking phenomena in domains with varying thickness is also of great im¬ 
portance. Ventricular fibrillation is a state of electrical anarchy in part of the 
heart that leads to rapid chaotic contractions of the heart. It is fatal unless 
a normal rhythm can be restored by defibrillation. One of the medical hy¬ 
potheses is that ventricular fibrillations are triggered by the generation of a 
circular excitation wave that are thus trapped in a region of the heart. These 
waves are not present in a healthy heart and an important medical problem 
is to understand how these waves can be triggered. Ashman and Hull m sug¬ 
gested that a myocardial infarct could create a circular area of cardiac tissue 
where excitation waves can propagate only in one way. Once an excitation 
wave would enter this area, it would be trapped and the propagation of the 
wave would trigger ventricular fibrillations. 

Hence an important question is to understand how excitation waves can 
propagate only in one direction in cardiac fibres. Cranefield pA] asserted that 
asymmetry of conduction in a normal cardiac fibre occurs due to variations 
of the diameter of the fibre. This asymmetry is strengthened by depression 
of excitability. 

Here too, mathematical models of excitation waves in the cardiac tissue 
may shed light on these phenomena. Cardiac cells as brain cells are excitable 
and the propagation of the depolarisation initialising the excitation wave 
can be modelled by a bistable reaction-diffusion equation as in (1.1). One of 
the important problems is thus to understand if the geometry of a fibre (or 
rather a fibre bundle) may trigger an asymmetry in the propagation of the 
travelling front and if geometrical properties of the domain may even block 
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propagation of travelling waves in one direction but not in the other. This 
is precisely what we establish here. Indeed, our findings show that in the 
presence of an abrupt transition, waves can be blocked in one direction (the 
direction in which the domain is widening) whereas they propagate in the 
opposite one (the direction in which it is narrowing). 

Related earlier works include numerical studies [23;, [3T] that supported 
this hypothesis. Grindrod and Lewis (2D studied the propagation of the 
excitation in a Purkinje fibre bundle using an eikonal equation obtained as 
an asymptotic limit of a model of propagation. They used biological values 
for all parameters and computed the approximated speed of the front. In the 
context of this asymptotic limit, they proved that both an abrupt increase 
in the diameter of the bundles as well as an heterogeneous depression of 
cardiac tissues may prevent propagation of fronts in one direction preserving 
the propagation in the other direction. 

Lastly equation (1.1) is also relevant to describe spatially structured pop¬ 
ulation dynamics. Bistable equations, of the type we study here, model 
populations that are subject to an Allee effect. One can think of the inva¬ 
sion by a population of fishes, through a straight into wide ocean. Invasions 
of species (plants, trees, animals) subject to Allee effect, going through an 
isthmus and then into a large area also lead to the kind of geometric effects 
we study here. 

In (SB], the authors studied numerically the evolution of a population in 
such a complex domain as it faces climate change. In this framework, the non¬ 
linearity in equation ( 1.1[ ) is space dependent: / = f(x , u ) with, say f(x, u ) = 
—7 u when \u\ > A (where 7 and A are positive constants). Furthermore, 
to model the effect of climate change one considers / = f(x — ct,u ) where 
c is the exogeneously given climate change velocity. The domain involves 
a transition region with an opening angle a so that when a is small, the 
transition is gradual whereas when a is close to 7t/2 the transition is abrupt. 
It is shown numerically in [38] that there exists a critical angle «o £ (0, 7 t/ 2) 
such that when a > ao there is extinction of the population, whereas if 
a < a 0 , the species survives. This property is analogous to what we prove 
here. Using our methods here we can derive the results when a is either cloe 
to 0 or to 7 t/ 2. Nevertheless, the complete description of [[38] with a specific 
critical angle is still an open problem. 
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Existence and uniqueness of the entire so¬ 
lution 


The proof of Theorem 1.4 follows directly from the work of the first author 


with Frangois Hamel and Hiroshi Matano on waves passing an obstacle |7]. 
We will give here the main points needed to adaptat the arguments but for 
a detailed proof we refer to [7, section 2.3-2.4-3]. The idea is to construct 


sub- and supersolutions of problem (1.1) for negative time t < T < 0 for 


some T < 0 and to use these sub- and supersolutions to construct the entire 
solution. For the uniqueness result one uses the monotonicity in time and a 
comparison principle. 


3.1 Construction of the entire solution 

The construction of the entire solution follows exactly the same steps as in 
section 2.3]. In our setting, the super- and subsolutions are defined by 


w + (t, x ) = 


and 


w ( t , x ) = 


2 </?(—ct — £(£)), if xi > 0, 

<p(xi — ct — £(£)) + v ?(—xi — ct — £(£)), if xi < 0, 


0, if X\ > 0, 

<p(x 1 - ct + £(£)) - <p(—x 1 - ct + £(£)), if Xi < 0, 


(3.1) 


(3.2) 


where the function £ is such that 

i(t) = Me x ^ ct+m \ £(—oo) = 0, 

with A = - (c + y/c 2 — 4/'(0)^ and M a positive constant to be determined. 
Thus, 


A ln (l -Mc-'e** 


for all t < In ( -f- 
~ Ac \M 


If we assume that ct + £(f) <0 which is true as soon as t < j~ c In j > if 
can be proved that there exist T < 0 and M > 0 large enough such that w + 


is a supersolution and w is a subsolution of problem (1.1) for all t <T and 
ieO. 
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Now the idea is to construct a sequence of solutions u n of the Cauchy 
problem defined for — n <t< +oo, such that u n converges toward an entire 
solution as n —> + 00 . Following [TJ section 2.4], we define u n to be the 
solution of (1.1) for all t > —n with the initial condition 

u n (—n,x ) = w~(—n,x). 

Since w~ < w + and using the comparison principle one proves that w~ < 
u n < w + and 

u n (t, x) > u n _i(t, x) for all t G [—n + 1, + 00 ), x G f l. 

Using the monotonicity of the sequence and parabolic estimates as n —)■ + 00 , 
we see that u n converges to an entire solution u of (1.1) for all t G R, x G 0, 
and 

x) < u(t, x) < w + (t, x) for all t G (— 00 , T\, x G U. 

Since ^(— 00 ) = 0, we have thus proved that there exists an entire solution 
to (1.1) such that 

| u(t, x) — ip(x 1 — ct ) | — > 0 as t —> —00 uniformly in iGfl. 


3.2 Time monotonicity 

As in [7], we will make use of the property that Ut > 0 for all t G 1 and 
x G U. This can be proved using the maximum principle and the fact that 
w~[ > 0. Hence 


(u n ) t (t,x) > 0 for all t G (—n, + 00 ), x G H. 


Letting n —> +00 and using the strong maximum principle, one concludes 
that 


u t (t, x) > 0 for all t G M, x G H. 


Since 0 < u < 1, using monotonicity and usual parabolic estimates, we obtain 
that u(t, .) —» Woo in Cf oc (Q) and is a solution of (1.10). 

In the following, we will need more informations on u t so for 0 < i) < |, 
we introduce 

Vt^t) = {1 G fl : rj < u(t, x) < 1 — r)}. 


19 






Lemma 3.1. For all rj G (0, |], there exists S > 0 such that 

u t (t,x ) > 5 for all t G (—oo,T v ] and x G 12, ? (t). (3.3) 

Proof of Lemma 3.1\ The proof of the lemma is the same as in [7] section 

3 except that the domain is changed. We will therefore give only the main 

ideas. 


First using (1.3), there exist T v G R and M n > 0 such that, 
Fl v (t) C {x G 12, \x\ — ct | < M v } C x cu for all 


oo < t < T v . 


(3.4) 


Now argue by contradiction and assume that there exist tk G (—oo, T v ] and 
Xk G 12 n (t) such that u t (tk,Xk ) —» 0 as k —> +oo. There are two possibilities: 

• If (tk)k is bounded, then Xk is also bounded and up to extraction of a 
subsequence, we can assume that tk —* t M and Xk —> x* so u t (t*, x*) = 0. 
But u t > 0 on 12 so x* G 912 and using the parabolic Hopf lemma, 
dvUt = ( d u u) t < 0 and we obtain a contradiction with the boundary 
condition. 


• If {tk)k is unbounded, then Xk is also unbounded which means here that 
xp*; tends to —oo. Defining Uk(t, x) = u{t + tk, x\ + x\ t k, x') and using 
parabolic estimates, we get that Uk —> u* as k —y +oo in C^(f xRxw). 
Then, we conclude as in the previous case. 

□ 


3.3 General comparison principle and uniqueness of 
the solution 


We now want to prove that the entire solution satisfying condition (1.3) is 
unique. In a sense, this result is an extension of the parabolic comparison 
principle to the case when the initial condition is given at t = — oo, which 
we now state in the following more general form. 


Lemma 3.2 (General comparison principle). Assume u and v are super- and 
subsolutions of (1.1) with limiting condition (1.3), that is: 


{ d t u — A u > f(u ) in R x 12, 
d u u >0 on M x 912, 
lim inf u(t,x) — ip(x — ct) > 0. 

t- 5—oo xefi 


and 


d t v — Av < f(v) in R x 12, 
d u v <0 on M x 912, 
lim sup x (t, x) — tp(x — ct) <0 


(3.5) 
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Then v < u for alltEM. and x G fl. 

This comparison principle indeed contains the uniqueness of the entire 
solution satisfying (1.3). 

Proof of Lemma 3.2 1 In the following proof, we will note Cw := w t — A w — 

/H- 

Choose rj > 0 small enough so that 

f'(s) < -7 for s G [- 277 , 2r/] U [1 - 277 ,1 + 277 ] (3.6) 

for some 7 > 0. For any e G (0, 77 ) we can find t £ G R such that 

u(t, x) — w(t, a;) < £ for — 00 < t < t £ and x G 

We let T* = min(T ?; — ae,t £ ). For any t 0 G (— 00 , T*], we define 

z(t, x ) = u(t + to + cre{l — e _7< ), x) + £e -7 * 

where the constant cr > 0 will be specified later. Then 

v(t 0 ,x) < z( 0,x) for x G Ll. (3.7) 

Let us prove that z is a supersolution in the time range t G [0, T* — f 0 ] if a 
is sufficiently large (independently of e). 

Cz > a r j£e^' yt u t — 7 £e _7f + f{u) — f{u + £e _7t ) 

= £e~’ yt (a r ju t — 7 — f(u + 9ee~ Jt )), 


for some 9 — 0(t, x) with 0 < 6 < 1 . 

If x G Tl v [t + t Q + ct£( 1 — e' 7 *)], then using the preceding lemma 3.1 
see that 


we 


Cz > ee ^fcT'yS — 7 


max f'(s )). 

0<s<l J V ” 


Therefore, Cz > 0 if a is chosen sufficiently large (independently of £ > 0). 
If x fL Cl v [t +1 0 + < 7£(1 — e -7 *)], we have u(t +1 0 + cr£(l — e^ 7 *), x) + Oee^ 1 * G 
[0, 277 ] U [1 — 77 ,1 + 77]. Consequently, f'{u + 9ee~ li: ) < 7 and 


Cz > ee 7 *(— 7 + 7 ) = 0. 


Applying the classical comparison principle for t G [0, T* — 1 0 ] and x G T2, we 
get 

v(to + t,x) < z(t, x) for t G [0, T* — t 0 \, x G ft. 
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Rewriting t 0 + t as t for simplicity of notation, we see that 

v(t, x) < u(t + ae(l — e -T'(^*°)) ; x ) + ee ~ 7 h-*°) f or t G [t 0l T*],x G fl 
where to is chosen arbitrarily in (—oo,T*]. Letting t Q —>■ —oo, we derive 

v(t, x) < u(t + ae, x) 

for all t G (— 00 , T*],x G ft and we can once again apply a classical compar¬ 
ison principle and conclude that this inequality holds for all 1 G R. Finally, 
letting e —> 0, we have proved v < u. □ 


3.4 More general domains 


The proofs we have sketched here actually work for more general domains. 
For instance we can consider bent cylindrical type domains in which we have 


existence and uniqueness of an entire solution of problem (1.1) satisfying 
(1.3). By slightly modifying the sub- and supersolutions w~ and w + for 
X\ < 0, we can extend these results to U-shaped cylindrical domains. Of 


course in these type of domains, condition (1.3) is satisfied for only one of 
the end of the cylinder near x\ —y — 00 . We also have the same existence 
and uniqueness proof for domains that only have one infinite “end” near 
X\ —» —00 and the remaining part is bounded. These different types of 
domains are illustrated in Figure |4} 





Figure 4: Examples of 2D domains for which existence and uniqueness of an 
entire solution of Problem (1.1), satisfying (1.3), hold. Left: domain with 
one end ; center, bent cylindrical type domain ; right: U-shaped cylindrical 
domain. 
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4 Propagation in domains with decreasing cross 
section 


In this section we consider propagation in the direction of a narrowing do¬ 
main. We formulate this condition by requiring the first component of the 
outward normal to be non-negative, see figure [5] for examples in 2D. Here we 


prove Theorem 1.6 






ii 




Figure 5: Examples of 2D domains concerned by theorem 


1.6 


We start by proving that the unique solution u of (1.1) such that (1.3) is 
satisfied cannot be blocked and propagates to 1 at large times. The proof is 
based on the following lemma. 


Lemma 4.1. Under the same assumption as in Theorem \1.6[ in particular 
that the domain D satisfies vfix) > 0 for all x 6 dTl, let (ip, c ) be the travelling 
wave solution defined by (1.4) and u be the unique solution of (1.1) such that 
(1.3) is satisfied. Then for allt^R and 

tp(x i — ct ) < u(t, x). 


Proof of Lemma (4.1). Define v(t,x) = <p(x\ — ct) for all ( t,x ) e M. x D, v 
satisfies the following problem: 


d t v — Av = f(v) in Ixll, 

d v v <0 on M x dQ, 

^v(t, x) — <p{xi — ct) 0 as t —> —oc uniformly in i6 0. 
Indeed, for all t G M, x G dtt, 

d u v(t,x) = <p'(x i — ct)ui(x). 


(4.1) 


23 



























As assumed in Theorem 1.6, v\{x) > 0 for all x G dil and since p is de¬ 
creasing, we see that d v v < 0 in M x <9hL Applying the general comparison 
principle of Lemma |3.2[ we get v < u. Q 


Next, we make precise the behaviour of the solution u at large time, i.e 
when t —» +cx). Since u is bounded between 0 and 1 and cp(x i — ct) —> 1 
when t — > —f-oo uniformly locally for x G fl, we have u(t,x) —> 1 uniformly 
locally for x G D. 

To complete the proof of Theorem 1.6, we now furthermore assume that 

\>N — (1 i w l ^ n , , r- iraV-l 


D fl {x G Xi > 1} = (l , +cxd) x cu r for a given / > 0, ay C 
prove that the ID speed c is the asymptotic spreading speed i.e. 


We will 


for all c > c 
for all c < c 


lim sup u(t, x) = 0 , 

t—>+0° Xl> Qt 

lim inf u(t,x) = 1 . 

£—>•+00 X\ <ct 


(4.2) 

(4.3) 


The lower bound for the speed of propagation in (4.3) immediately 
from u(t,x ) > ip(x i — ct) for all ( t,x ) GRxh 

We now turn to the proof of ( 4.2[ ), which gives an upper bound 
speed of propagation. Let us fix c > c. For all r] G (0, |), introduce 
Figure[ 6 ]) a smooth function in [ 77 ,1+?/] such that /(?/) = f (9) = /(I + 
/ < 0 on ( 77 ,0 ), f > 0 on ( 0 , 1 + 77 ) and 


results 

for the 
fn (see 
77 ) = 0 , 


fv = f in [ 2 ? 7,1 - 77 ], 
fn>f in [ 77 , 277 ) U (1,1+ 77 ]. 


(4.4) 


Consider the 1-D travelling front solution (p v of speed c n associated with f rn 
that is ((p v , c v ) satisfies the equation 


1 ++) + + fr,M = 0 for all zGR, 

[+??( 00 ) — 1 T 77 , ^(+ 00 ) = 77 . 

It is well know that c v —> c as 77 — >■ 0. We denote 770 > 0 the threshold such 
that c v < c for all 77 G (0, 770 ). Let 77 G (0, 770 ) be fixed. Since | u(t, x) — <p(x\ — 
ct) | —> 0 as t —> —00 and lim +0O cp — 0, there exist T < 0 and L > 0 such 
that for all t < T and x\ > L, u(t, x) < 77 . 
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Figure 6: Example of function /„ that satisfies (4.4). 


Let us consider the following problem 


'v t (t, x ) — A v(t, x ) = f(v(t, x)), for t G [T, +oo[, x G [l, +oo[xuy, 

d u v(t, x) = 0, for t G [T, +oo[, x G [/, +oo[x<9o; r , 

v(t, l, x') = u(t, l, x'), for t G [T, +oo[, x' G 0 J r 

k n(T, x) — u(T, x), for x G [/,+oo[xo; r . 

(4.6) 

It is constructed in such a way that u is a solution of it and v = tp v (xi — c v t—£,) 
is a supersolution if £ G K is chosen large enough so that — £) > 1. 
According to the comparison principle, v > u. Thus 


lim sup u(t,x ) < lim sup v(t,x) = 7] 

t-++oo xi> g t t^+OO X1> £t 


since c > c v . As tj can be chosen arbitrarily small in ( 0 , 770 ), we infer from 
this that 

lim supw(f,x) = 0 

7-H-oo Xl> £ t 


and (4.2) is proved. 

Remark 4.2. The same kind of result can be obtained for smooth pertur¬ 
bations of straight cylinders. Precisely, using techniques from ra , one can 
easily prove that if is a family of domains that converges in the C 2,a 
topology to a straight cylinder as £ —» 0, then u £ the solution of (1.1) with 
condition (1.3) converges to 1 everywhere. Moreover if v\(x) < 0 for all 
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x G dfl and fl £ D {xi > 1} = (l, +oo) x u j £)V , the solution u £ propagates with 
the 1-D speed c. 


5 Blocking by a narrow passage 


In this section we prove Theorem |1.8 

There are earlier results on the existence of non constant stable steady 
states for bistable reaction-diffusion equations with Neumann boundary con¬ 
ditions when there is a narrow passage in a number of cases. Matano [[29] 
first showed the existence of non-constant stable solutions in the case of a 
bounded domain having the shape of an hourglass. The first author of the 
present paper together with Hamel and Matano [7] established an analogous 
result for an exterior domain having a narrow passage to a confined region. 
The paper [7] exhibits a non constant stable steady state in an exterior do¬ 
main with narrow passage where the area in which the solution is close to 0 
is bounded whereas the area where the solution is close to 1 is unbounded. 

Our proof is inspired from |7J. Here, we construct a non constant stable 
steady state for domains where the two areas are unbounded (the one where 
the solution is close to 0 and the one where the solution is close to 1). Stability 
is obtained since the solution is constructed by minimization of an energy. 

Actually, we derive a more precise result here. We will show that once 
the measure of fl D {b < aq < b + 1} is fixed, there exists £ > 0 small enough 
such that if |fl D {a < aq < 6}| < £ then there exists a supersolution w of 
the elliptic problem on fl D {aq > a} such that w — 1 on O fl {aq = a} and 
w goes to 0 in fl, as aq —* +oo. 

Thus, we consider the reduced problem in fl' = {x G fl, aq > a}, (see 
Figure [F]). 


Aw + f(w) = 0 in fY, 
d u w = 0 on d£l'\{xi = a}, 

w = 1 on {x\ = a}, 


(5.1) 


Such a w blocks the propagation of the solution u of (1.1) satisfying (1.3). In¬ 


deed, by the comparison principle, u(t, x) < w(x) in fl', whence u a0 (x\, x') —>■ 
0 as X\ —y Too. 


To find a solution of (5.1) we proceed in two steps. We first restrict our 


analysis to a bounded domain Q' R = fl' fl {aq < R}, R > b + 1 (see figure 
[7]) and prove in section 5.1 that there exists a function 0 < Wr < 1 in f 1] 


R 
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Figure 7: Illustrations of the various domains used in the proof of Theorem 1.8 for 
two different 12 satisfying assumptions ,12[ ) and ( |1.11[ ) in Theorem |1. 8 1 12 is the initial 
domain, 12' = {x G 12, X\ > a}, fl' R = {x £ 12, a < X\ < i?}, 12 t = {a: £ 12, a < X\ < b} 
and 12 fi = (cc £ 12, b < X\ < R} 


solution of: 


'Aw + f(w) 
d u w = 0 

< 

w = 1 
= 0 


0 in 0' fi , 

on T = dtt' R \({xi = a} U {xi = R}), 
on {xi = a}, 
on {or = R}. 


(5.2) 


Then we prove in section |5.2| that w R —> as R —> Too with —> 0 

as X\ —> Too and conclude, using the general comparison principle, that 
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the solution u of (1.1)-(1.3) is blocked by in the right hand part of the 
domain. 


5.1 Reduced problem 


For D C M. N , we introduce the energy functional defined for all w € H 1 (D), 
by 

JdH = [ hvw\ 2 + F(w)dx, (5.3) 

where F(t) = J 1 f(s)ds. Note that using (1.6) and extending / linearly 
outside [0,1], F(t) > 0 for all tel and that F has the shape illustrated in 
Figure [8j 



Figure 8: The function F(s) = J g f(t)dt, where / is a bistable function for s € [0,1] and 
is extended linearly outside [0,1]. 


Let Fl = {x e dFl' R , x\ = a} and = (i 6 9Fl' R , x\ = R}, we define 

Hi o = jw e il 1 (n / R ), w E 1 on V L and w = 0 on j. 

Finding a local minimizer of Jqi r in an open subset of H\ 0 yields a stable 
solution of (5.2). 

We start by proving that the constant function wq = 0 on = {x € 
12, b < x i < R} is a local minimum of in Ft 1 (12#). 








Proposition 5.1. There exist 5 > 0 and a > 0 such that for all R > b + 1, 
for all w E with ||rc||#i(n fl ) < S one has 


Ja„(w) > Jn K (tOo)+a||TO||gi (0ll) . 


(5.4) 


Proof of Proposition \5J\ We first observe that, 
Jn R (w) = 

JQr 


^|Vw | 2 + F(0) + F'(0)w + ^^w 2 + rj(w)w 2 , (5.5) 


where r/ : R —>■ R is a continuous function such that r](s ) —)■ 0 as s —> 0. This 
can also be written as 


Jn R (w) = Jn R (wo) + [ \\Vw? - + ri(w)w 2 . 

■P-h> z z 


(5.6) 


To exploit (5.6), we require the following lemma: 


Lemma 5.2. If rj : R —>■ R is the function defined in (5.5), we have 


r](z)z“ 


'^R 


< h{ z )\\ z \\ 2 Hi(n R )> 


(5.7) 


where n(z) ->■ 0 as ||^||i?i(n K ) ->• 0 


Proof of Lemma 5. By extending / linearly at infinity, we may assume that 
F is quadratic at infinity, which implies that rj is a bounded function from 
R to R, i.e there exists C G R such that \rj(s) | < C for all s € R. Let us 
fix p > 0 sufficiently small so that f) C L p+2 (Hr) independently of 

R > b + 1. Since rj is bounded and rj(s) —* 0 as s — * 0, for all 7 > 0 there 
exists £ > 0 such that 

\v(t)\ < 1 + & P . (5.8) 

Hence, 


r](z)z z 


'Of 


— 7lkll/fi(n H ) + ^11^11^+2(0^)- 


(5.9) 


By using Sobolev embedding we infer that 



— (7 + 


II 2 

\\m(n R )- 


(5.10) 


Since 7 can be chosen as small as we want, this inequality yields Lemma 5.2 


□ 
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As /'(0) < 0, we can fix a = \ min(|, —Qp-) > 0. Thus, 




J|Vw| 2 - > 2a\\w\\ 2 H i iaR) 


and from the previous lemma, there exists 6 > 0 small enough such that for 


all — 8, one has r](w ) < a. Substituting in (5.6) yields 

Jn R (w) > Jn R (w 0 ) + Oi\\w\\ 2 m ^ R) . (5.11) 

and proposition 5T is proved. □ 

Now going back to the domain Q' R , we extend the function Wq by 

f° r Xl ^ l a ib), X G 12, 


W 0 (x) = Wo(x-i) ■ = 


(5.12) 


0 for x\ G [5, R ], x G 12. 

The function wq is in H\ 0 . The next step is to prove the following proposition: 

Proposition 5.3. There exists 6 > 0, such that for all w G H\ G with ||tc — 
woll/rpoy = 8 one has 

Jn'(w) > Jndw 0 ). (5.13) 


Proof of Proposition \ 5.3\ In the following, all constants will be denoted gener- 
ically by C. For all w G H\ 0 , 

Jn'Jw) - J Q > (w o) = Jn T {w) - Jn T (w 0 ) + Jn R (w) - Jn R (w 0 ). 


(i) 


( 2 ) 


If ||w-Wo||ffi(ny = <2, one has that ||w-Wo||/ri(n fl ) = IMItf 1 ^*) - 5 and one 
can use Proposition 5.1 to prove that the second term (2) above is bounded 
from below, i.e 

Jn R (w) - Jn R (w 0 ) > a\\w - w 0 \\ 2 h i^ r) . 

For term (1), as /(0) = 0, /'(0) < 0 and F(l) = 0, there exists k > 0 such 
that 


F(s) > n(s — l) 2 , 


30 












for all s6l (see Figure [8]). This implies that 

Jd(w) > v\\w- 1\\ 2 h i { d) , (5.14) 

for all D compact domain in tt' R , w G i/ 1 (Zl), where v = min{AC, |}. Moreover 

Jn T ( w o) = / ^|Vw 0 | 2 + F(w 0 ) < C\Q T \ and |0 T | < e. (5.15) 
Jn T 4 

and ( |5l5| ) , we obtain 

Jqt(w ) - Jn T {w 0 ) > u\\w - l||^i (nr) - Ce. 

Then we also use 

u \\ w — l||/p(f2 T ) ^ ~j\\ w — u ’o||f/ i (n T ) — v\\wq — l\\ 2 H i { a T ) 

> ^\\w - w 0 \\ 2 h i^ t) - Ce. 

Therefore there exists 7 > 0 such that 

Jttjw) ~ Jn' R {w 0) > 7||w - w 0 ||^i (n y - Ce. 

Taking £ > 0 small enough, we get 


From (5.14 


J, 


n' R {w) - Jn' R (w 0 ) > 0 . 


This completes the proof of Proposition 5.3 


□ 


From Proposition |5.3[ we get the existence of a local minimizers wr of the 
energy functional J^ R that belongs to //,] ] such that ||wr — u>o|| < S. This 
yields a stable solution of (5.2). For S > 0 small enough, we have 0 < wr < 1 
using the maximum principle. 


Remark 5.4. Note that in this proof we could actually choose any fonction 
g(x') on ff D {07 = a } and make the same construction and prove that there 
exists a solution of 


'Aw + f(w) = 0, in Q' r , 

d u w = 0 on T = dQ' R \({xi = a} U {x± = R}), 

w = g(x') on {07 = a}, 

= 0 on {07 = R}. 

for £ small enough. 

The next step will be to study the behaviour as R —> +00. Let us recall 
that S is fixed for all R> b + 1. 
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5.2 Construction of a particular supersolution 


In this section we construct a particular supersolution of (1.1). 
with the following Proposition: 


We start 


Proposition 5.5. Let wr be the minimizer of the energy functional Jq/ r 
defined in the previous section. When R +oo, up to a subsequence, wr 
converges to a solution Woo of 


{ A Woo + f{w oo) 

(dpW oo 0 

Woe = 1 


0 in fl', 

on T = <9fl'\{xi = a}, 
on {xi = a}, 


(5.16) 


such that Woo —)• 0 as x\ —* -Poo. 


Proof of Proposition 5. As 0 < wr < 1 and using Schauder estimates, 
there exists R n 
and Woo 

Woo —> 0 as xi —> -poo. Using Fatou’s Lemma: 


—> -poo as n —> -poo such that wr u — > Woo in Cf oc as n —* +oo 
is a solution of (5.16). It remains to prove that the limit Woo satisfies 


IK - Woo|K (n0 < liminf||K - w Rn )t { _ 1 <Xl< R n} \\ 2 m{iil) < S 2 


Then arguing by contradiction, we assume that there exist rj > 0 and a 
sequence (x n = (:r n)1 , :r(J) neN such that x n G hi for all n G FT, x n% i — > -poo 
as n — > -poo and Woo{x n ) > r/, for all n G N. As Woo € Cf„ c , it implies that 
|Vwoo| < C on every compact set K. For all x G B(x n , ffe) fl hi, 

|Woo (x) - Woo (x n ) I < max || Vwoo|| \x - x n \ 
x€B(xn,fy)ntt 

So Woo{x) > |, for all x G B(x n , fl fl, for all n G N. This yields 

11 Woo - Wo\\l 2 {K) > 2 n fl| x ( number of x n G AT). 

As fl is assumed to be uniformly C 1 , ||woo — WolK#) > S for K large enough, 
which is impossible. EJ 

We now extend w^ by 1 outside fl', i.e 


Woo (x) 


Woo(x) if x G fl', 

1 if x G fl\fl'. 


(5.17) 
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Then is a supersolution of the parabolic problem 


v t — Av = f(v) in fi, 
d v v = 0 on dd, 


(5.18) 


Moreover lim inf w oc — (f(xi — ct) > 0. Indeed for x\ < a, w 00 — tp(xi — ct) = 

t—>■—oo x£fl 

1 — tp(xi — ct) > 0 for all t G M. and for X\ > a, lim ip(xi — ct) = 0 

£->-—oc 


uniformly in x\ > a. Using the generalised maximum principle 3.2[ we 
obtain u(t,x ) < Woo(^) for all tel and iGfl. As u is increasing in t and 
bounded between 0 and Woo^x), u(t,x ) — > u^x ) as f —> +oo for all a; G U 
and passing to the limit in time, u^x) < w^x) for all x G fh We have thus 


proved Theorem 1. 


6 Partial propagation in domains containing 
a sufficiently large cylinder 

In this section we investigate the propagation properties of u when U contains 
an infinite cylinder in the xi-clirection with large enough radius. We will 
prove Theorem 1.10 We assume here that M x B' R C U where R > Ro for 


some R 0 which we now define. 

We consider the following semi-linear elliptic equation with Dirichlet con¬ 
dition in the ball Br of radius R in with N > 2: 


-Az = f(z ) in B r , 


z = 0 


on dB 


(6.1) 


R- 


For large R, positive solutions can be found by minimization of the energy 
functional: 

«,) := /' / 

JBr 2 


f(s)ds 1 dx 


in Hq(Br). It is known [8] that there exists a threshold Rq > 0 such that 


for R < Rq, there is no positive solution of (|6.1|) while for all R > Rq, there 


exists at least one positive solution of (6.1). Moreover, for all R > Rq, there 


exists a maximum positive solution which we denote by z = zr in the sequel. 
Incidentally, from a result of Rabinowitz EH. it can be inferred that for all 


R > Rq, Problem (6.1) has at least two distinct solutions. 
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The maximum solution satisfies 0 < Zr < 1 in Br and is radially sym¬ 
metric and decreasing with respect to \x\ [THj. We extend Zr by 0 outside 
Br and denote Zr(x) = wr(\x\) where wr is defined on M + . In particular, 
we know that sup zr = 2 r( 0 ) = wr(0) > 6. 

Using the function wr we will now derive the following estimate from 
below on for all R> R 0 : 

Lemma 6.1. For all x = (. xi,x ') G 12, we have 

Uocix) >w R (\x'\). ( 6 . 2 ) 

Proof. This is obtained by a simple sliding argument. We know that for fixed 
R, we have max vjr < 1. Therefore, for large enough a > 0, we have u^xi — 
a,x') > maxvjR > iur(|x|) for all x = (x\ , x r ) G B R , hence everywhere in 12. 
Let u a (xi,x') := u oc (xi — a,x'). Both functions u a and zr are solutions of 
the same elliptic equation in Br and Zr — 0 while u a > 0 on the boundary of 
Br. We also know that u a > zr for large enough a. We now start decreasing 
a. From the strong maximum principle we infer that the function zr has to 
stay below u a for all a G M. Consequently, 

Uoo(xi,x') > maxzfl(xi + a, x') = Wr(\x'\) for all x± G M, x' G 

agR 

which proves the claim. □ 

Next, we prove that u 00 is bounded from below by a positive number up 
to the boundary of 1 x B' R . This is stated in the following lemma. 

Lemma 6.2. There exists 6 > 0 such that 

u 00 (xi,x') > wr(\x'\) + 5 forx£M.xB' R . (6.3) 

Proof. Assume by contradiction that for all n G FT* there exists x n = (aq in , x' n ) G 
M x B' r such that 

Uoc(Xl,m x'n) < Wr(\x'J) + -. 

Consider v n (x ) = Moo(^i,n + %i, %') ~ w^d^l) a function defined on M x B' R . 
Then, up to extraction, we can assume that x' n —> x' 0 , v n —> v in C Rc (M. x B' R ) 
and v verifies 

v > 0, u( 0 ,a: / 0 ) = 0 , v ^ 0 , 

—Av = c(x)v with c bounded, 
v > 0 in dB’ R . 
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The last positivity property is obtained from u 00 > 0 in O and by using that 
Moo satisfies the Neuman boundary condition, and passing to the limit in it, 
if need be. This contradicts the strong maximum principle. Q 


As a conclusion, inf > 6 and the proof of Theorem 


IxB' 


complete. 


1.10 


is thereby 

□ 


Remark 6.3. Note that, Lemma 6.2 actually proves a somewhat stronger 
result than the one stated in Theorem 11.101 


In order to conclude that the state u = 1 invades the whole domain, one 
has to require further assumptions on the geometry of fl . Indeed in domains 
as illustrated in Figure [9] there is propagation in the sense given above but the 
propagation could be blocked when entering some parts of the domain. One 
can find in [7] or [L0] arguments leading to partial propagation in exterior 
domains and examples of domains where the propagation is actually only 
partial. 



Xi 


Figure 9: Example of domain in dimension 2 where propagation could be blocked in the 
small area on the top of 17. 


We now turn to this question of finding conditions under which the inva¬ 
sion by 1 is complete. Since u(t, ■) —> Uoo locally uniformly in as t —> oo, 
the limit is a solution of the following stationary problem 

{ -Aiioo = /(Moo) in D,, 

d u u 0 c = 0 on dft, (6.4) 

Uoo(x) —* 1 as xi —* —oo, 
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We are thus led to finding geometric assumptions on that imply that 

Uoo — 1 • 

It is interesting to note that this is related to a classical result of Matano 
[29] and Casten and Holland [12] regarding the following non-linear elliptic 
Neumann boundary value problems in bounded domains D: 


—A u = f{u) in D, 
d u u = 0 on dD. 


(6.5) 


These authors have shown that if D is a convex and bounded smooth domain, 
then the only stable solutions are constants. It is a natural question to ask 
whether there are extensions of this result to unbounded domains. What 
further conditions imposed on the solution and which geometrical conditions 
on the unbounded domain D garantee that the only stable solutions are 
constants? Our results on complete propagation are partial answers to this 
question for the solution Uoo- More general formulations are an interesting 
open problem. 

We have shown in Theorem 11.81 that there exists a stable non constant 
solution of the previous stationary problem for domain containing a transition 
area of small enough measure. We have also shown in Theorem L6 that 
Uoo = 1 for propagation in the direction of a decreasing domain. In the next 
sections we derive two different types of geometrical conditions that garantee 
the complete invasion by 1. 


7 Complete propagation for domains that are 
‘star-shaped” with respect to an axis 




The object of this section is to prove the result of complete invasion of Theo¬ 
rem 1.12 that is Uoo = 1 in fh We consider a domain H such that fl D lx B', 


_ 

with R > Rq, Rq > 0 defined in section [6} Moreover we assume that at each 
point on the boundary x = {x\,x') G dfl, the outward unit normal v makes 
a non-negative angle with the direction x'. To prove the result, we introduce 
here a new method which rests on sliding cylinders. 

As in section [6j for all R > Rq, we denote z = zr the maximal positive 
solution of 

-Az = f(z) in B Rl 
z = 0 on BBr. 
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It satisfies 0 < z R < 1 in B R and is radially symmetric and decreasing with 
respect to |x|. We extend z R by 0 outside B R . We denote z R (x) = iur(|x|) 
where w R is defined on M + . In particular, we know that sup z R = zr(0) = 
w R { 0 ) > 9. 

By shifting the axis of the cylinder in orthogonal directions and rotating 
the cylinder, we will now construct from w R a new axisymmetric function 
that is a subsolution. Let h > 0 a given real and e G a given unit vector 
defining a direction in M^ -1 . We define, for all x' G l^ -1 : 


i>h,e{ x ') := w r(\ x ' ~ M) 


so that 'i^h, e ( x ') > 0 on B R (he) := B R + he, and 

ip h (x') := max ^ he (x'). 

e £§ N ~ 2 


(7.1) 


Note that, since tc^(|a; / |) is axisymmetric and decreasing, iph(x') = w R (\x'— 
h Jx'l I) f° r x ' 0- Hence cis also axisymmetric but not decreasing away from 
the origin anymore. Moreover iph(x') > 0 if and only if h — R < \x'\ < h + R. 
Finally since iph is a supremum of a family of subsolutions, it is a generalised 
subsolution, as defined in [6]. But we will only use classical subsolutions. 
We consider it to be defined on all of M. N , so that it is independent of the 
variable x\. 

By Lemma 6.2[ we know that u 00 (x ) > w R (\x'\) = c po( x> ) for all x G fl. 
We claim that u 0Q (x) > ifih( x ') f° r all x G O and all h > 0. Hence u^x) > 
ip\ x >\{x') = w _ r ( 0 ) > 9 and we will prove later on that this implies Uoo = 1 . 

So let us assume by contradiction that there exist h 0 > 0 and a;o G H 
such that Moo(^o) < ^(^o)- We define 


h* = sup{/i > 0 , \/h' G [0, h] Vx G H Uoc(x) > (7.2) 


Clearly, h* > 0 is well defined and finite. By continuity, we know that 
Uoo{ x ) > y>h*{y) f° r all igL 

By definition of h*, there exists a sequence h n \ h* and a sequence of 
point x n G H such that u^Xn) < <fih n ( x ' n ) where x n = (xi in ,x' n ). Without 
loss of generality, we can assume that the sequence is bounded. Indeed, if 
the sequence (x n ) n is not bounded, then we can shift the origin by (x n ) n so 
that now Xi >n = 0 but the domain and the function change with n. We can 
pass to the limit on these by a standard compactness argument. 
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Therefore, we can assume that the sequence (x n ) n converges to some point 
p = (0 ,p') G In the limit, we get Moo(p) < ph*{p') and since Uoo > ifh * 

(p r ) = < fih*{p')■ First, we know from Lemma 6.2, that 
(xi, 0) > max p + S for all X\ 


everywhere, we get u, 
there exists 5 > 0 such that u. 


G M, so that it is 

impossible that p' = 0. It is also impossible that p G Indeed, there exists 
some e G S> N ~ 2 such that u^p') = Ph*(p') = V J h*,e(p r ) > 0 and u 00 > 'iph”,e 
everywhere. In view of the strong maximum principle (recall that u a0 — 
satisfies some linear elliptic equation, since u 00 and i^h*,e satisfy the same 
semi-linear equation), we infer that u 00 = But this is impossible as 

Vv, e vanishes somewhere. 

It remains to consider the case that p G dfl. Since u a0 > 0, we know that 
h* — R < \p'\ < h* + R. Let us show that, necessarily, h* < \p'\ < h* + R. 
Indeed, suppose to the contrary that h* — R < \p'\ < h*. Recall that p' ^ 0. 
In the region h* — R < r < h*, the function (fh is decreasing with h. Thus, 
for small enough h* — h > 0 we see that Ph(p') > <Ph*(p') — u oo(p)- But 
this is in contradiction with the definition of h* that requires > <ph for 
such an h. We have reached a contradiction and, therefore, we know that 
h* < \p'\ < h* + R. 

We have now reached a situation where u 00 > ph* everywhere, Uoo(p) = 
= wr(\p' - ^*]^[|) = ™ R (\\P'\ ~ h*\). 

We compute the outward normal derivative using the fact that Wr is 
decreasing to get 


d v <p h *(jf) = - W R {\p' 




The term jR comes from the derivative of \p' — h* with respect to p' and 
the information that h* < \p'\ < h* + R. The last inequality follows from our 
assumption 


v 


p 

W- > 0 


|p'I 


on 


on 


The proof is thereby complete. 

Now we want to prove that = 1 in Assuming that ^ 1, there 
exists a sequence (x n ) n in such that u OQ (x n ) —* infu^x) G [iw^(0),l) as 

xef i 

n —> Too. We recall that Tr(0) > 6. Letting u n (x) := Uoo(x n T x) and 
D n = {x G M. N , x T x n G 0}, we know that u n u and D n —> D as 
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n —» +cx) such that 


—A u = /(w) in D, 
d u u = 0 on dD. 


If 0 G dD, it contradicts Hopf’s Lemma and if 0 G D, then Att(O) < 0, but 
f(u( 0)) > 0 which is impossible. D 


Remarks 7.1. 


The assumption (1.13) is a star-shaped type property. It 


essentially says that any section in planes orthogonal to the X\- axis is 
star shaped with respect to the trace of the axis in this plane. Hence 
the name “star shaped with respect to an axis” 


• We can also treat the case when ends to the right”. Namely H is a 
straight cylinder as x\ —» — oo, and say is bounded to the right. 




Figure 10: Examples of domains that satisfy the sliding cylinder assumption in dimension 
2 . 



£ 



Figure 11: Examples of domains that do not satisfy the sliding cylinder assumption in 
dimension 2. 


This proposition proves the complete propagation of 1 in a large variety 
of domains but examples of domains that are not covered by this theorem 
are presented in figure 11 Some domains, as the figure on the left of Figure 
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11, that do not satisfy starshape assumption with respect to an axis never¬ 


theless satisfy the conditions of Theorem |1.13| and thus satisfy the complete 
propagation property as we will now prove in the next section. 


8 Complete propagation in the direction of 
an increasing domain 


In this section we prove Theorem 1.13 
following assumptions: 


Thus we assume that satisfies the 


x B' r CO 


(8.1a) 


where R > Ri and R\ is defined below. We further assume that there exists 
L > 0 such that , 



:= i( x h x ') G fi, X\ > L} is convex, 

(8.1b) 

and that there exists C > 0 such that 



{x G ff, X\ < L + i?} C (— oo, L + R) x Be , 

(8.1c) 

and 

z/i(x) < 0 for all x G dVt with x\ < L + R. 

(8. Id) 


We will first prove that Uoo is close to 1 in 

r L+R O fi {(xi, x') G M x xi> L + R } 

using a sliding method. Then we will extend this estimate to all of using 
a backward travelling front as subsolution and finally we will conclude that 

Uoo = 1 * 

As in section [6j we denote z = Zr the maximal positive solution of 

<-Az = f(z) in B a , 

I z = 0 on 8Br, 


We recall that it satisfies 0 < Zr < 1 in Br, and is radially symmetric and 
decreasing with respect to \x\. We extend zr by 0 outside Br and denote 
Zr(x) = iur(|o;|) where Wr is dehned on M + . In particular, we know that 
supyjj = Zr(0) = Wr(0) > 6. We further recall from Lemma 6.1 that for all 
x = (xi, x') G fl, we have u^x) > wr(\x'\). 
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Figure 12: Domains 12 that satisfy (8.1) 


By a result of [8] we know that this maximal solution wr converges to 1 
locally uniformly as R +oo. Thus, there exist R\ > Rq, where Rq > 0 is 
defined in section [6j and S > 0 such that wr( 0) — 2 5 > f3 for all R > Ri, 
with (3 > 9 defined by (1.7) in the Introduction. Let rj > 0 be such that 
wr(|o;|) > wr(0) — 5 for all x G B v . We will prove that > wr(0) — 5 in 
^l+r using the following definition. 

Definition 8.1 (Sliding ball assumption). We say that x G 12 verifies the 
sliding ball assumption if there exists a G B v (x ) 


x° G 


and 


7 : [ 0 , 1 ] —» Q continuous, such that 7 ( 0 ) = (x°, 0 , ..., 0 ), 7 ( 1 ) = a, 

and V£g[0, 1) V b G Br^^)) D (b — 7 (f)) • v > 0 (H2) 

where v is the outward unit normal at b G dQ. 


Assumption (H2) means that there exists a continuous path 7 from a a 


point near x to (x°, 0 ,..., 0 ) a point of the xi-axis, such that for all x t on this 
path B R (x t ) fl 12 is star shaped with respect to the center x t . 
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Lemma 8.2. For all x Eft that satisfies the sliding ball assumption, we have 

(x) > w R (0) - 5. 


Un 


Proof. The proof is the same as the one of Lemma 6.1 with the novelty here 
that we slide the maximal solution w R along the path between the xi-axis 
and a and compare it to the solution u 00 . 

So we consider x E hi that satisfies the sliding ball assumption and an 
associated path 7 from (x°,0) on the xi-axis to a E B. q (x) satisfying the 
properties in the definition. Let w R be the maximal solution of the Dirichlet 
boundary value problem (8.2) in B R (x). We know from Lemma 6.1 that 


w 


(x0,0) 
R 


< Uoq, so we can define 

t* = sup{f G [ 0 , 1 ], Vf G [0 ,t\ W ^ < Moo}- 

We will prove by contradiction that t* = 1. Otherwise we define v = u c 
uPr ^ and once again we have 

—An = c(x)v in with c bounded, v > 0 , v ^ 0 , 


there exists x G B R ( 7 ( 2 *)) such that v(x) = 0. 


By the strong maximum principle, it is impossible that x G B R ( 7 ( 2 *)) flO so 
x G B R (-y(t*)) fl dfl and from Hopf’s Lemma, we have d v v(x) < 0. But 

x — 7 (£*) 


d u v(x) = -d„(w R (\x - 7 (t*)|)) = \w' R (\x - 7 (t* 


• V 


\x- 7(01 

since w R is decreasing. But this contradicts the definition of the sliding ball 
assumption. 

We thus get that u 0Q > w R and since x G B v (a), we have Moo(^) > 
w R (0 ) — 5 by definition of r/ and S. □ 

We claim that if x — fx\,x') G O verifies x\ > L + R, then 

Moo(:r) > w R (0) - 5. 

Indeed let us choose x E fl with x\ > L + R. Thanks to the previous lemma, 
it is sufficient to prove that x satisfies the sliding ball assumption. Since fl r L 
is convex, as assumed in (8.1), we define a = x, xq = x\ and 7 (t) = (xi, tx') G 
fl r L for all t G [0,1]. Assume by contradiction that there exist t G [0,1) and 
b G B r ( 7(t)) D dfl such that (b — 7 (t)) ■ v < 0. Then b and 7 (t) G fl r L which 
is convex so for any £ G [0,1], (1 — e)b + £7 (t) = b — e(b — 7 (f)) G fl r L but 
this contradicts the definition of the outward unit normal, for £ small. 
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Remark 8.3. Note that > wr(0) — 5>9 in kV L+R implies that u oc (x) —>■ 1 
as X\ —* Too. 

Now we will prove that Moo(x) > iur(O) — 28 > 9 for all x G hi and thus 
Moo = 1 in hi following the same proof than in the previous section. This 
ends the proof of Theorem |1.13 

Let us define f$ G C' 1 , 1 ([0, Wr(0) — 25}) as follows: 


fs 0 ) 


f(s) Vs G [0, /?], 

0 for s = w_r( 0 ) — 25, 


(8.3) 


with 0 < fs < / in (0, ?xr( 0) — 2<5). And consider (<^, 5 , c^) the travelling wave 
solution between 0 and Wr( 0 ) — 25, invading 0 from the right, i.e (ip 5 , c$) is 
solution of 


-<$ - C M = fs(<Ps) hi R, 

ps(~ 00) = 0, ^5(Too) = w R (0) - 25, c s < 0. 


We have u^x) —>■ 1 as x± — * —00 so we can fix — M < 0 such that Uoo(x) > 
Wr( 0) — 5 on kl l _ M := (x G hi, Xi < —M}. Now since u 0 0 > Wr( 0) — 8 > 
wr( 0) — 25 > ps also on {x G hi, x± > L + R}, we have > e on hi for e 
sufficiently small and thus we can translate p$ on the right such that, 


Vx G hi, u 00 (x) > p s (x 1 ). 

We can now apply a parabolic comparison principle to and v(t,x) = 
p&(x 1 — ct ) on M + x {x G hi, — M < x\ < L + R}. Indeed / is locally 
Lipschitz and 

on x {(xi, x') G hi, —M < x± < L + R} , 

on x {(xi, x') G hi, —M < x± < L + R} , 

for t > 0 and x G <9hl with — M < x\ < L T R, 

for t > 0 and x G hi with Xi = L T R or Xi = —M, 

for x G hi with — M < x\ < L T R. 

This yields Woo(x) > p(x 1 — cst ) for all t > 0 and x G hi with — M < x± < 

L T R and passing to the limit for t —* 00 , we obtain u 00 (x) > wr(0) — 2<5 on 
{x G hi, — M < X\ < L T R} and thus on all hi. 

Finally we conclude that = 1 as in the previous section. 

Remark 8.4. This proof can be adapted for more general domains than the 
ones considered in Theorem 1.13. The clue is to prove that the sliding ball 


dt^oo AMqo f (u 0 o) 

d t v s - A v s < f(v s ) 
d^u^t, x) = 0 > d v v(t, x) 
Uoo(t,x) > v(t,x), 
uoo( 0 ,x) > v( 0 ,x), 


assumption can be verified for any point of the domain far to the right. 
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